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. Abstract. Statistical invariance of Wiener increments under SO(n) rotations 

^ I provides a notion of gauge transformation of state-dependent Brownian motion. We 

■ show that the stochastic dynamics of non gauge-invariant systems is not unambiguously 

I , defined. They typically do not relax to equilibrium steady states even in the absence 

' of extenal forces. Assuming both coordinate covariance and gauge invariance, we 

. derive a second-order Langevin equation with state-dependent diffusion matrix and 

vanishing environmental forces. It differs from previous proposals but nevertheless 
entails the Einstein relation, a Maxwellian conditional steady state for the velocities, 
and the equipartition theorem. The over-damping limit leads to a stochastic differential 
' equation in state space that cannot be interpreted as a pure differential (Ito, 

^ ■ Stratonovich or else). At odds with the latter interpretations, the corresponding 

Fokker-Planck equation admits an equilibrium steady state; a detailed comparison 
with other theories of state-dependent diffusion is carried out. We propose this as a 
. theory of diffusion in a heat bath with varying temperature. Besides equilibrium, a 

^ I crucial experimental signature is the non-uniform steady spatial distribution. 
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1. Introduction 



Nonequilibrium thermodynamics is often modeled by means of stochastic differential 
^ ■ equations (SDE's), describing possible trajectories in an open system's state space, 
and by diffusion equations for the evolution of ensembles [H [21 [3] . While the paradigm 
started by enhancing dynamical equations with an external source of noise, independent 
of the system's state, soon there arose the necessity of considering situations where the 
intensity and direction of the environmental disturbance depends on the state of the 
system [U El Ej. From a different perspective, mathematicians had extensively studied 
SDE's on manifolds, where a state-dependent diffusion coefficient is often associated to 
a Riemannian metric (see Refs. [ZIIH] for recent monographs). At the merge between the 
two approaches Van Kampen argued that "The question of the existence and correct 
form of equations describing Brownian motion on a manifold cannot be answered by 
mathematics alone, but requires a study of the underlying physics" [9]. 
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1.1. Scope and motivation 

While statistical properties of Brownian motion are expected to only depend on the 
diffusion matrix and the physical drift, loosely speaking the "square root" of the diffusion 
matrix appears in a SDE. In one dimension this is unambiguous, while in n dimensions 
there are many different ways to take such square root, differing one from another by 
local internal rotations. We will require the physical description to be independent 
of this choice. When this doesn't occur, for a given diffusion matrix one can write 
down SDE's that have very different dynamics and thermodynamics. In particular most 
choices lead to nonequilibrium steady states even in the absence of external forces: An 
unobservable statistical property oddly becomes a nonequilibrium drive. 

Taking up Van Kampen's invitation, in this paper we make a proposal for 
state-dependent diffusion equations based on the physical symmetry principle of 
General Covariance, by which we mean: (1) Covariance under external coordinate 
transformations; (2) Invariance under local internal rotations (gauge invariance). The 
covariance of nonequilibrium statistical mechanics under coordinate transformations has 
been extensively studied by Graham [10] (also Ref.[TT] for applications to semiclassical 
methods in field theory). Departing from Graham's analysis, we further assume 
that statistical properties of the theory are invariant under a class of internal gauge 
transformations under rotations of the noise components, taking inspiration from the 
frame field approach to General Relativity. However, we do not derive relativistic 
equations for diffusion in strong gravitational fields. 

Recently we witness an increasing interest in geometric methods in nonequilibrium 
statistical mechanics [121 [131 EH CSl [HI [H] , with emphasis often given to the equilibrium 
vs. nonequilibrium character of steady states. This paradigm will allow us to better 
appreciate the physical nature of our construction. It is our conviction that the theory 
is relevant to diffusion in environments with nonuniform temperature, as we briefly 
anticipate in the next paragraph and argue more thoroughly in a parallel publication 



1.2. Models and results 

In this paper we consider stochastic models of systems subject state-dependent diffusion, 
viscous damping and no external forces, unless otherwise stated. We do not consider 
velocity-depedent diffusion, as is done in Ref . [6] . To this purpose second order Langevin- 
type equations serve the purpose better than do first order SDE's. In fact, the latter 
are plagued with the problem of interpretation, which makes "vanishing forces" a fuzzy 
concept. Moreover, second order equations go onstage on a geometrically neat and tidy 
tangent bundle. 

In one dimension, the following equation with state-dependent diffusion coefficient 
and vanishing environmental forces is frequently encountered [5l [191 120] 




(1) 
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where y/2Ct is Gaussian white noise, a{x) = ae{xY is a state-dependent diffusion 
coefficient and 7(x) is state-dependent viscosity. We will later refer to this theory 
as "state-dependent viscosity". Assuming the state-dependent Einstein relation 7(x) = 
/3(t(x), the theory leads to a Maxwellian steady state P*(x,w) oc exp— /3f^/2, with 
uniform spatial density. Experimental evidence by Langon |21], Volpe [221 EH] and 
coworkers seems to indicate that Eq.([T]) describes the motion of an overdamped colloidal 
particle when the bath has uniform temperature, while displaying a nonuniform density 
due to hydrodynamic interactions with confining walls. 

Less studied [Ml §1.3.2.1, §4.1.2.2] [251 |26l |27] is the following equation 

xt = -jxt + V2ae{xt)Ct- (2) 

It was originally proposed by Van Kampen [2B] to realize a theory settling Landauer's 
physical intuitions [29] about systems with varying temperature, and therefore we will 
later refer to this theory as "state-dependent temperature". In fact, it is obtained from 
the state-independent Langevin equation by simply replacing /3 with f3{x) = e{xY'j/a. 
It differs from Eq.([T]) in one respect: The viscosity 7 is not itself state-dependent. It leads 
to a nonuniform steady distribution of particles p*{x) oc P{x), in accordance with the 
phenomenon of thermophoresis which drives particles along the temperature gradient, 
from hotter to colder zones. 

We advance yet another equation: 

Xt = -jxt + [d^ loge(xt)] x^ + V2ae{xt)Ct. (3) 

As the latter theory, it involves a state-independent viscosity. However, an additional 
quadratic term in the velocity appears — somewhat analogous to Rayleigh's drag 
force. The theory accomodates Einstein's relation and a Maxwellian steady conditional 
probability 7r*(f |x) = P* {x , v) / p* (x) for velocities at a given position, while the spatial 
density is generally not uniform, dispensing with the equal-a-priori postulate. In Ref. [18] 
we argue that this might be a valid alternative theory for systems with nonuniform 
temperature, even more adherent to Landauer's thermodynamic picture. In the present 
paper we discuss the geometrical properties of the theory and of its n- dimensional 
generalization. 

In the overdamping limit 7, cr — )• 00 at fixed (3, the viscous damping term is 
responsible for the reduction of Langevin-type equations to first order equations over 
state-space, where the inertia x of the particle can be neglected. Notoriously, the latter 
are subject to antique [30l |31] and ongoing [5l[20l[2ll[22l[23l[32l|33l|34| debates 
on the "correct" interpretation of the stochastic differential: Ito, Fisk-Stratonovich, 
Hanggi-Klimontovich, or else? Since different stochastic differentials are mapped one 
into another by adding effective drift terms, the problem of interpretation is equivalent 
to the choice of drift in a given fixed convention; throughout this paper we choose to 
work with Ito calculus. We talk of pure SDE when the equation has no drift in some 
convention. Remarkably, the generally covariant first-order SDE stemming from our 
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theory is not pure, as it requires and effective drift in all conventions, but for the one- 
dimensional case, where it reduces to the Stratonovich differential. It rather coincides 
with the mathematicians' definition of Brownian motion on a manifold. 

One remarkable property of our construction is that the steady state of the 
generally covariant reduced Fokker-Planck equation always satisfies detailed balance, 
with vanishing steady currents, while other approaches might lead to nonequilibrium 
steady states even in the absence of environmental forces, with steady currents driven 
by inhomogeneity of the diffusion coefficient. In particular this is always the case for all 
pure stochastic differentials. 



1.3. Plan and notation 

The paper is organized as follows. In Sec. [2] we discuss the statistical meaning of gauge 
invariance under rotations and derive the generally covariant Langevin equation, the 
corresponding Kramers generator and its steady state. Sec. [3] is devoted to the Fokker- 
Planck equation and to the corresponding first order SDE deriving from the overdamping 



limit, which is carried out in Appendix B (for a slightly more general theory). We also 
define covariant thermodynamic forces, currents, entropy and entropy production. In 
Sec. m we compare our results with other theories of state-dependent diffusion and 



provide comparative tables in Appendix D Finally we draw conclusions in Sec. [5l 



In the following, x, x' are local coordinates on an n-dimensional manifold M, which 
is assumed to be compact without boundary to avoid problems with conservation of 
probability; coordinates of a stochastic trajectory; Explicit dependences on 

coordinates and time are often omitted; Indices a,b,c and i,j,k label respectively 
coordinate and euclidean directions; The Einstein convention on index contraction is 
used; Spatial derivatives are abbreviated with da = d/dx"", Tensor indices are raised 
and lowered with the metric, Xa = QabX^, internal indices with the euclidean metric 
Xi = 5ijXK Fields and probabilities are assumed to be smooth, so that there is no 
nonequilibrium driving due to topological defects such as sinks and sources. 



2. Brownian motion 



2.1. Gauge invariance under rotations 

We work with n independent sources Q of Brownian white noise. They are the 
formal time derivatives of n independent Wiener processes {Wl)t- Wiener increments 
dWl = Qdt have null average, are uncorrelated at different times and in different internal 
directions, 

(dWi) = 0, {dWi dWi,) = 6'^6{t - t')dt, (4) 
and are normally distributed with covariance 6^^dt: 

Prob {dWD oc exp — ^— -J -. (5) 
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All three properties are invariant under rotations of the Wiener increments, dW^ — )■ 
R^jdW\ where Rij a special orthogonal matrix, 

R'^R'Ak = Sji, (6) 

with unit determinant R = det{Rij)ij = 1. The reader should pay attention to the 
collocation of indices to avoid confusion between a matrix and its transpose/inverse. 
For example, RjR\ = SI while R^jR\ = {R^)\- Invariance under global rotations is 
obvious for Eq.Q and is easily proven for Eqs.Q from the linearity of the average. 

The gauge principle suggests to consider the analogous local symmetry. The 
intuitive picture is as follows. Wiener increments are attached to each point of state 
space. As a Brownian particle passes by xt at time t it receives a kick. While single 
trajectories do change if we rotate the Wiener increment at that point by a rotation 
matrix R^j{xt), the statistics of kicks remains unchanged. In fact, linear combinations 
of independent multivariate normal variables have multivariate normal distribution, with 
average 

{RYxt)dW^) = {R'^{xt)){dWi) = 0, (7) 

and covariance matrix 

{R\{xt)dWt^ R^i{xt)dWl) = {R\{xt)R^^{xt)){dWt^dWl) = 6'^6{t - t')dt. (8) 

In both cases we used the properties of Ito calculus, which we conventionally work 
with. More precisely, any product h{xt)dW^ of a state function by a stochastic 
increment requires an interpretation about when the state function is evaluated within 
the time interval. The interpretation stems from the underlying discretization procedure. 
The trajectory is sampled at discrete times (tK)^, when it occupies states (xk)k- At 
intermediate times t G [t^j^K+i) one approximates h according to h{xt) = h{x^ + 
— Xk)), where: a = (Ito) has a causalistic, non-anticipatory character; a = 1/2 
(Fisk-Stratonovich) is middle-point and preserves the rules of calculus; a = 1 (Hanggi- 
Klimontovich) has a finalistic character. We emphasize that the choice of Ito calculus is 
purely conventional and does not affect the generality of the results. More importantly, 
this choice is not related to the dispute that we hinted at in the introduction. 

The advantage of Itb calculus is that the statistical properties of Wiener increments 
are independent of the state Xt where they are nested, hence we are allowed to 
factor as we did in Eqs. fl7|8|l . What these equations then tell us is that the rotated 
Wiener increments have the same statistical properties as the starting ones: They 
are statistically indistinguishable. This is a gauge transformation in that it acts on 
the internal properties of the system — we do not measure the stochastic increments 
themselves — and it depends different iably on the state. 

Since Wiener increments enjoy this symmetry, we expect that any theory whose 
statistical properties are dictated by the statistics of Wiener increments also enjoys 
this symmetry. As we will point out, just like theories that are not coordinate 
covariant depend on the particular coordinates chosen, theories which do not satisfy 
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gauge symmetry have very different properties according to how the symmetry is made 
to brealc. Gauge invariance of overdamped Brownian motion has also been recently 
discussed in Ref.[35]. 



2.2. Generally covariant Langevin equation 

While the the manifold of states M needs not be physical space, for sake of visualization 
the Langevin equation 



is thought to describe the motion of a Brownian particle with mass m. The trajectory 
is driven by an effective drift F", is damped by a viscous force —70;°, and receives 
stochastic kicks from the Q's. The frame mixes n independent Wiener increments to 
yield n coordinate increments. It determines intensity and direction of random kicks, in 
a way that depends on the state visited by the trajectory. Correlations are characterized 
by the diffusion matrix 

g'^'ix) = enx)e%x)6'^. (10) 

We assume that (e^)j^a is an invertible matrix, that is, the n frame vectors = (e")a 
are linearly independent. Its inverse (e^)i,a is called the coframe. It follows that 
the diffusion matrix is positive definite, with inverse {gab)a,b and inverse determinant 
g = det{gab)a,b > 0. 

We now perform an invertible coordinate transformation x — )■ x'{x) and, by the 
principle of covariance, we assume that the trajectory x[ = x'{xt) obeys an equation in 
the same form as Eq.Q: 



x"" =v\ v" = --^v" + + V2ae^ C- (H) 

In theory, we should resort to the rules of Ito calculus to evaluate differentials; in 
practice, for this kind of equation stochastic and standard calculus coincide, and there 
is no issue concerning the interpretation of the stochastic differential. From the first 
transformed equation x""' = v°''daX°'' we read off the transformation law for the velocity 
y"-' = v°'daX"'' , which behaves like a vector sitting in the tangent space at x. We then 
take the stochastic differential of v°'\ using the (Ito-)Leibniz rule 

dv''' = v^'ddax'"' + dv^dax""' = [dad^.x'"' v^v^ + dax'''F'')dt + y/2dax'''e1Cdt, (12) 

which leads to ^ 

= -T^el = ^ + v^v \ ^ , . 13 

* dx'' dx" dx'^dx^ 

The frame's components, at fixed i, transform like vectors, while the inverse 

diffusion matrix is twice contravariant, gab = ga'b'daX'^' d^x^' . Therefore, the latter bears 

the properties of a metric, with invariant volume element ^/gdx. Since gabe-le^^ = 6ij, the 

frame vectors form a basis of orthonormal vectors for the tangent space at x, whence 
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the interpretation as orthonormal reference frame, point-by-point over the manifold. 
In general, such reference frames are anholonomic (noncoordinate), meaning that the 
matrix ef is not the Jacobian of a coordinate transformation that makes = 6^ and 
gab the flat euclidean metric, but at one point on the manifold at a time. The frame 
makes its appearance in General Relativity, where it is known as vielbein (counting 
dimensions in german, viel = zwei, drei, vier, etc.). For example, the vierbein maps the 
local observer's 4D fiat Minkowsky metric — three sticks and a clock, following Einstein's 
physical intuition — into a semi-Riemmanian metric on the spacetime manifold. 

As is well known, the inhomogeneity developed by the drift term in Eq.f llSp qualifies 

it as 

= -r>V, (14) 

where F^^ are the Christoffel symbols of the Levi-Civita connection, i.e. the unique 
torsionless connection compatible with metric Qab, 

2TI = g'^" {dbgdc + d^g^b - d^gcb) ■ (15) 

However, this is not the only viable choice, as 

= v^'eld^e'i (16) 

produces the same inhomogeneous term. As we will see later, in the overdamping 
limit the latter leads to Graham's theory while the former leads to the mathematicians' 
Brownian motion on a manifold. 

For a sensible choice between the two, we invoke the principle of gauge invariance. 
We assume that the theory is invariant under internal SO(n)-transformations of the 
reference frames, acting on the i,j, k indices 

etix) ^ -etix) = e^^{x)R\{x), e\{x) ^ e^x) = e^.(x)i?;(x), (17) 

where the transormation law for the coframe follows from e"e^ = 5^. Notice that the 
metric is invariant under local coordinate transformations, cjab = gab, and so are the 
Christoffel coefficients. In General Relativity, gauge transformations have the physical 
meaning of rotations of the observer's orthogonal reference frame. As we argued in the 
dedicated paragraph, Q and = R\{x)Q are statistically indistinguishable, hence the 
stochastic term is invariant under e^Q ^tQ — ^iQ- 

Given that the stochastic term is statistically preserved by gauge transformations, 
it remains to check that the deterministic term is. As a rule of thumb, it should be 
expressible only in terms of the metric and of its derivatives. This is the case for Fi, 
while it is not for F2. In fact, after a gauge transformation: 

6F^ = 0, 6F^ = eteiv^v'R/d,R''i. (18) 

Therefore we choose Fi, obtaining 



v'' + F>V = --fv^ + V2^e^C- 



(19) 
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In the deterministic case 7 = a = one recognizes the geodesic equation. Eq. (fT^ is 
our proposal for a Langevin equation with state-dependent diffusion matrix. On a flat 
manifold, in coordinates where the metric is euclidean and the Christoffel coefficients 
vanish, it returns the ordinary Langevin equation. 

At the time of submitting, the author came to knowledge that Eq. (fT9l) can already 
be found in a work by Kleinert and Shabanov [36], who discussed its generalization 
to connections with torsion. Very recently it has also drawn the attention of Castro- 
Villareal and coworkers [37], who derived it using an extrinsic, rather than intrinsic, 
approach. 



2.3. Kramers equation and its steady state 

While the n-dimensional process {xt)t is not a Markovian diffusive process, the 2n- 
dimensional process Zt = {xt, ik) on the tangent bundle TM is. By a standard 
procedure, the Kramers equation is obtained by writing down the diffusion equation 
for the probability density Pt{z) and expressing it in in terms of {x, v): 

dP 8 

The steady state is given by 

P*{x,v) = 3Mew-2-g,,{x)v%\ (21) 

It can be checked by direct substitution; the calculation entails an interesting balancing 
of Christoffel terms T abcV"'v'^v'^ and 2r^^w'^ from v to x derivatives. 

Since the metric has physical dimensions [a;]~^, the viscosity and the diffusion 
constant have dimensions respectively [7] = [t\~^ and [a] = [t]~^. We define an energy, 

H{x,v) = ^gUx)v'^v\ (22) 

where we introduced the Compton wavelength of the Brownian particle Ac = 2'Kh/mc 
for sake of dimensional consistency. In Appendix A it is shown that the average energy 
Et = {H{xt,Vt)) is a Lyapunov functional, exponentially decreasing to the steady value 
nl3~'^/2, where 

7/(T = f3{27ihf/mc'^. (23) 

The identification of /3 with an inverse invariant temperature is suggested by the 
equipartition theorem at the Gibbs steady state P* oc exp —PH, yielding the Einstein 
relation. Upon coordinate transformations, the energy is a scalar, the temperature is 
invariant and the probability density transforms like a squared volume density (like g), 
consistently with the fact that Pt{x, v)dx A dv ought to be an invariant volume element 
to preserve normalization. 

Finally, we calculate the steady spatial density by integrating over velocities: 



hoo , 



p*{x)= I P*{x,v)dv= ' — (24) 

iMy9{x')dx' 
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where = (mc^/27r/3/i^)"/^ J^^ g{x')dx' . Brownian particles tend to occupy space 
according to the volume element. More probable states occur where noise is less intense. 
In fact, if we modulate the intensity of noise without modifying the angles between frame 
vectors, by performing a rescaling Cj — )■ Ajej (no index contraction), the determinant 
of the metric scales like g — )■ g/Wi^i- This is physically consistent, as more intense 
noise kicks particles further away, making permanence at a state briefer. Modifying the 
angles, g increases when the frame vectors are made less and less linearly independent, 
shooting to infinity when (e^)i,a becomes degenerate. 

The conditional probability density 7r*(t>|x) = P*{x,v)/p*{x), describing how 
velocities in the tangent space T^M distribute at a fixed x, is Maxwellian. However, since 
X and V terms in P* do not factor, the velocity probability density 11* (t^) = jj^jP*dx 
will generally not be Gaussian. This is a distinctive signature of the theory. 



3. Overdamped dynamics 

3.1. Fokker-Planck equation 

The overdamping (Smoluchovski) limit of the Langevin equation is achieved by sending 
7, a — !■ oo at fixed /3, in which case the damping effect of the viscous force is so strong 
that the inertial term drops. The limit must be taken with great care; details of a 



rigorous derivation are reported in Appendix B, where the Kramers equation is shown 



to reduce to the following Fokker-Planck equation 

Pt = da [g""' {d,vt - Vt d, log ^)\ (25) 
for the spatial probability density 

Pt{x) = / dvPt{x,v). (26) 



Introducing the covariant derivative Vq associated to the Levi-Civita connection, defined 
by VflX'' = daX'' + r'^^X'^, Eq.( l25l) can be recast in manifestly covariant form 

p = Amp - V, (/» , (27) 

where pt = Pt/ y/g is a scalar and = VqV" is the Laplace-Beltrami operator. For sake 
of generality we momentarily accounted for an environmental force vector f"". In deriving 
Eq.( l27|) . we used the fact that the covariant derivative of a scalar coincides with its 
partial derivative and that the covariant gradient is given by WaX'^ = y/g ^^da{y/g X""). 
When = 0, Eq. (l271) returns the well-known definition of Brownian motion on a 
manifold in the specialized mathematics literature (up to a factor 1/2) [71 Ch.3] [H 
§4.2.2]; It was also considered by Batrouni et al. in Ref. [38] with particular interest to 
the stochastic quantization of lattice gauge theories, and by Castro Villarreal [39], who 
worked out the small-time/low-curvature dynamics. 
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3.2. Smoluchovski equation 

The SDE that corresponds to Eq.( l25l) is given by 

i'^ = r + /eff + V2e?C, (28) 

where again we included the environmental force vector, and the effective drift reads 

= -g"^ri = 8,9^" + log ^. (29) 

Eq.f l28|) can also be derived from first principles in much the same way as the generally 
covariant Langevin equation, keeping into account that by the rules of Ito calculus the 
effective drift transforms like 

f:;, = d,x'^'{f:s+g'^%d,x^'). (so) 

In short, one performs a Taylor expansion of x'{xt) to second second order in dxt, plugs 
Eq. (1251) in and retains terms of order dt, keeping in mind that squared Wiener increments 
(dW^)^ equal dt with certainty. 

Interestingly, when = 0, for no value of a does Eq. fl28l) coincide with a pure SDE 

Xt = V2e'}{xt + adxt)Ct, (31) 

as there always appears an effective drift term, but in one dimension, where it reduces 
to the Stratonovich differential. We notice in passing that the Ito differential leads 
to a gauge invariant, but not coordinate-covariant SDE, while the converse occurs 
for the Stratonovich convention. In all cases, upon either a coordinate or a gauge 
transformation, such differentials usually gain an additional effective drift. More 
interestingly, pure SDE's in Eq.f lST]) usually lead to nonequilibrium steady states, as 
we sparsely argue in Sec. HI 

3.3. Geometric insights 

The promotion of global symmetries to local ones by virtue of the gauge principle usually 
entails the appearance of a new field, which compensates the transformation law and 
introduces a new interaction. The gauge field for S0(r2) transformations is the spin 
connection Uaij, which generates an infinitesimal rotation as one displaces vectors along 
a given direction, viz. {ujaij)i,j is an element in the Lie algebra so(n) of skew-symmetric 
nxn matrices. Here we show how our theory writes nicely in terms of the spin connection 
and of the closely related anholonomicity vector. This is the only truly geometrical part 
of the paper; the reader might safely skip it. Nevertheless, we maintain the general 
inaccurate tone depicting geometry as index gymnastics. 

Since fields are tensors, we should hereby resort to Stratonovich calculus, which 
preserves the rules of differential calculus. Eq.( l28|) . with = 0, now reads 



(32) 
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where h{xt) o Q = h{xt + dxt/2)Q is the Stratonovich differential. In the first passage 
we added a drift term to switch convention. After some manipulation one obtains the 
second equation, where 

u^a = e^idael - d,ei) (33) 

is called the anholonomicity (co)vector. After a gauge transformation the 
anholonomicity vector gains an inhomogeneous term which on average balances exactly 
the one developed by the Stratonovich differential 

6 {uja) = - (e^BidjaC) = S {V2et{x,) o q'j . (34) 

Notice that the anholonomicity vector has the typical transformation law G~^da G of 
gauge fields, where G is a group element. Therefore, the anholonomicity vector is 
the gauge field of generally covariant stochastic thermodynamics. It will be present 
in all expressions where rotation invariance requires adjustment. In a way, the 
anholonomicity vector introduces an effective interaction between the Brownian particle 
and the environment which bends its motion. 

Let's see how the theory writes in terms of the spin connection. The latter is 
used to displace objects with internal indices along the manifold rigidly, according 
to VaX* = <9aX* + ul^jX\ while spatial indices are displaced with the Levi-Civita 
connection. A condition for a metric-preseving connection is that the covariant 
derivative of the coframe shall vanish, 

Vacl = daei - Tl^el + ul^ei = 0. (35) 

Antisymmetrizing, and recalling that the Levi-Civita connection is torsionless, = 0, 
we can express the anholonomicity vector as 

Ua = ujl^eie\, (36) 

where we used a;*j = S^^Uaij = 0. Finally we can write Eq.( l32l) in a particularly simple 
form in terms of the spin connection and of the the velocity vl = e^(xf) o x'}: 

v' + ey; = V2C. (37) 

Let us briefly discuss the geometric significance of the anholonomicity vector. It 
measures to what extent the orthogonal basis fails to be a coordinate basis. If there 
exists a collection of potentials such that = daU^, then u"" vanishes. Performing 
coordinate transformation x"" — )■ u^{x), the transformed metric is easily seen to be 
euclidean. In this case, the manifold is fiat, and it has vanishing Riemann curvature. In 
general, the Riemann curvature of a manifold does not vanish, so that for such systems 
it is impossible to find coordinates where the anholonomicity vector vanishes. We skip 
the case where the manifold is only locally flat. We only observe that (local) flatness 
has no thermodynamic relevance to our theory, as far as we can appreciate. 
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3.4- Covariant entropy and entropy production 

Entropy and entropy production are central to macroscopic thermostatics and 
thermodynamics. We give covariant analogues in first-order theoy, while a complete 
thermodynamic treatment is deferred to a future publication. Covariant equilibrium 
thermostatics has also been recently discussed in Ref . |35] . 

From an information-theoretic point of view, the uniform distribution maximizes 
the Gibbs-Shannn entropy, a measure of ignorance with respect to a reference prior, given 
by the uniform distribution itself. The Gibbs-Shannon entropy notoriously develops an 
inhomogenous term when one keeps the uniform prior fixed through a coordinate change, 
6S = — J plogdetdx' dx. It can be avoided by transforming the prior as well, which 
will therefore not be uniform in the new coordinates. 

This to say that it does not make much sense to pick the uniform prior in some 
preferred coordinates in the first instance. Our theory comes with its own preferred 
prior, namely the volume element y^. A natural covariant candidate as a measure of 
excess information is then 

= j dx y/g plogp, (38) 

where p = pj ^ is a scalar. The above expression is relative entropy with respect to ^Jg. 
Though not uniform (a matter of concern for some authors [IQ]), p* is analogous to the 
microcanonical ensemble, with equal volumes in state space having equal probabilities. 
There is an obvious circularity in this argument, since the choice of ^/gdx as the 
preferred volume element is subjective, as is equiprobability. In fact, y/g bears more than 
a resemblance with Jeffrey's prior in Bayesian inference, whose objective vs. subjective 
character is still debated [H]. 

The author recently argued that subjective priors have full citizenship in NESM 
[T7] . A change of prior does modify the thermostatics of the system but not its 
thermodynamics, which is fully encoded in the concept of entropy production, which 
we now introduce. Let us light up external forces for sake of generality. We introduce 
covariant current and affinities, 

r = rp - V V, Aa = fa- Vfe log p, (39) 

and the entropy production 

(T = j dx^rAa = j p{dx) gabA'^A' = (A^) . (40) 

The equalities on the r.h.s. show that the entropy production is non-negative. It 
vanishes at the steady state when the condition of detailed balance holds, that is when 
there exists a potential such that fa = — Va0, with steady state p = y^exp— 0. 
Notice that our free theory, with fa = 0, satisfies detailed balance. We mention that a 
covariant necessary condition for detailed balance is 

Vafb - Vbfa = dafb - d^fa = 0. (41) 
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Even if this condition is satisfied though, the force might not be globally a gradient, 
igniting topological steady currents. Notice that the covariant condition is equivalent 
to the noncovariant condition (second passage) by virtue of the torsionlessness of 
the connection. We are tempted to speculate that torsion might therefore be a 
nonequilibrium drive, but we leave this to future inquiry. 

The entropy production can be expressed in terms of the time derivative of the 
relative entropy, which quantifies the internal entropy production, and an additional 
environmental term, a = —D ^ + {A°'fa), where 

(A^fa) = I p{dx) r {fa + da log ^ - \ogp) . (42) 

The second term between curved parenthesis can be shifted from the environmental 
to the internal term, returning the standard Gibbs-Shannon entropy and the standard 
environmental term usually considered in the literature [3]. The covariant splitting 
of the entropy production is one of many possibilities that are equivalent up to the 
choice of reference prior, which is an additional fundamental symmetry of nonequilibrium 
thermodynamics, deeply mingled with coordinate covariance. 



4. Comparison of models 

We consider several models of state-dependent diffusion with vanishing exernal forces 
that can be found in the literature, and compare them with (A) our model: (B) Graham's 
covariant theory of nonequilibrium statistical mechanics, (C) Diffusion in nonuniform 
temperature, (D) Diffusion with state-dependent viscosity. In the overdamping limit (B) 
and (C) lead respectively to Stratonovich and Ito pure SDE's. For sake of completeness, 
we also consider (E) the pure SDE with posticipatory character (Klimontovich 
convention). Finally we include two more general theories: (F) a noncovariant theory 
with state-dependent viscosity and temperature, and (G) its covariant counterpart, 
a theory where the diffusion matrix and the metric need not necessarily coincide. 
The latter includes both our theory and diffusion with state-dependent viscosity as 



special cases. In Appendix B and Appendix C the overdamping limit for the latter 



two theories are performed. All others follow as special cases. We resume all these 
theories' properties in tables IDll and ID2I and through the interdependence diagram in 



Appendix D Notice that (G) preceeds (D) and (F) in the discussion. In this section we 
will set (2TihY /mc^ = 1. 



4-1- (B) Second order Graham theory 

Suppose to replace Fi by F2 in our theory. By taking the overdamping limit, along the 
same lines exposed in Appendix B one would obtain Graham's first-order theory [TO] . 
Therefore, we can already conclude that the theory is coordinate covariant but it is not 
gauge invariant. The derivation leads to the following Stratonovich pure differential 
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equation 

i« = e'l o C = b'^e%e] + ^2 e^C ■ (43) 
The corresponding manifestly covariant FP equation is 

P = V,(VV-a;», (44) 

where we remind that bj°- is the anholonomicity vector, defined in Eq. fl33p . As 
expected, the theory is not gauge invariant since the anholonomicity vector is not. 
Hence, knowledge of the covariance matrix is not sufficient for determining the 
theory. All statistical properties depend on the choice of frames including the 
equilibrium/no nequilibrium character. If for some choice of frame the theory satisfies 
detailed balance, uJa = — Va0, since the inhomogeneous term in Eq.f l34p is not generally 
a gradient, the gauge transformed theory might not satisfy detailed balance. As an 
example, consider a two-dimensional theory with = 5" and = in some gauge 
and coordinates. In two dimensions Rl{x) = Rl{d{x)) is a rotation matrix by an angle 
^{x). To first order in § the transformed anholonomicity vector (71,72) = {—d2'&,di'd) 
is certainly not a gradient. 

This problem with Graham's theory was also recognized in Ref. |42] . 



4-2. (C) State- dependent temperature 

To our knowledge, theories with state-dependent diffusion coefficient and uniform 
viscosity have been considered in the literature only in one dimension [251 ISHl |271 128] , 
where they have been used as models of systems with varying temperature. Let us 
generalize to n dimensions: 



In this quite general case temperature might not be isotropic, besides being 
inhomogeneous: It depends on the direction of motion of the Brownian particle. Given 
the considerations developed in this paper, the theory is gauge invariant but it is not 
coordinate covariant. The corresponding Kramers generator 

does not admit a locally Maxwellian steady state, as one can appreciate by plugging 
e'x.Y>[—hah{x)v°'v^/2 + k{x)] as an ansatz. First and third order terms in the velocities 
can only vanish if /c = and daQbc = 0. The exact form of the steady state of Eq.( l46i) is 
not easily computable to our knowledge. 

As for Graham's theory, in general the steady state does not satisfy detailed balance. 



The nontrivial overdamping limit is performed in detail in Appendix C leading to the 
Fokker-Planck equation in Ito form 

P = -daf = dadt [g'^'p] . (47) 
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The steady current is j* = —p*{gaddc9^'^ + da\ogp*). For it to vanish, hence for the 
steady state to be equihbrium, it is necessary that QaddcO'^'^ be an exact differential. 
However, this is not generally the case. Detailed balance does hold when temperature 
is isotropic, by which we mean that there exist coordinates such that 

gabix) oc T{xy^Sab, (48) 

In this case the metric is said to be conformally flat. 

In a parallel publication we argue on physical grounds that our model should be 
considered as a valid alternative to diffusion in nonuniform temperature [18]. 



4-3. (G) Metric /diffusion theory 

One peculiarity of our theory is that the inverse diffusion matrix coincides with the 
metric: Noise determines the form of deterministic geodesic motion. This provides an 
elegant theory, but it is not strictly mandatory for a generally covariant theory. 

We might introduce a new nondegenerate metric hab-, with Christoffel coefficients 
VL^^. Let the diffusion matrix g"-^ be a twice covariant symmetric tensor, possibly different 
from the inverse metric h"-^ . The equation 

x"" + Vtl^ x^x" = --f^x^ + V2ae'^C\ (49) 

with g"-'' = Sijclj^el, is generally covariant when 7^ is a tensor. We assume the Einstein 
relation 7^ = Pag°''^hcb- The diffusion matrix and the metric govern respectively the 
noisy and the deterministic behavior. The viscosity coefficient is a tensorial object, 
once covariant and once contravariant. It is a hybrid of the diffusion matrix and of 
the metric, interpolating between stochastic and deterministic behavior. The Kramers 
equation affords the equilibrium steady state 

P* = h exp -habv''v^/2 (50) 

Notice that the metric, and not the diffusion matrix, determines the form of the steady 
state. Our theory is recovered by setting hab = gab- The overdamping limit is performed 
in 



Appendix B , yielding 

P = da [h'''gcdh'''{dbp-pdtlogVh)] . (51) 

4-4- (D) State- dependent viscosity 

Theories with state-dependent diffusion and viscosity are widely studied. See Ref. [S] for 
a review. The n-dimensional case, with a rigorous treatment of the overdamping limit, 
is discussed in Ref. [19] , while Ref. [13] gives a microscopic derivation in the spirit of the 
Caldeira-Leggett model and Refs. discuss the path integral formulation. 

The free model is a special case of the above tensor viscosity model when hab = ^ab 
and fi^^ = 0. This identification can only hold in one given set of coordinates, hence the 
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theory is not covariant. Although, it is gauge invariant. The steady state is the usual 
Maxwell-Boltzmann state, displaying no dependence on the spatial variable: In a way, 
it is blind to the inhomogeneity of the medium. In the overdamping limit, the theory 
leads to the SDE 

xt = d,g''' + V2eKi. (52) 

The Fokker-Planck generator is p = da [g^'^dbp) , sometimes named after Fourier and 
Fick. The steady state always satisfies detailed balance. In one dimension Eg. (152!) 
reduces to a pure SDE in Klimontovich interpretation. 



4-5. (E) Pure posticipatory SDE 

Theories (B) and (C) led respectively to the Stratonovich and to the Ito differential 
in the overdamping limit. For sake of completeness, we also consider the pure SDE in 
Klimontovich form: 

= V2e'^{xt + dxt)C = 2 5'^e%e] + V2e1{xt)C. (53) 

The theory is neither covariant nor gauge invariant. Moreover, the drift term arising 
in the overdamping limit, as for the Stratonovich and the Ito case, is not generally a 
gradient, so the theory does not accomodate detailed balance as an essential feature. 
It is remarkable that, in higher dimensions, none of the pure stochastic differentials 
generally affords equihbrium steady states. 



4-6. (F) State- dependent viscosity and temperature. 

To complete the collection, we recall that Van Kampen had considered a 1-dimensional 
generalization of state-dependent temperature and state-dependent viscosity to the case 
where both are state dependent [16]. Jayannavar and Mahato [17] claimed to have 
derived the theory from a microscopic model of state-dependent interactions with a 
bath, although it appears that the state-dependence of the temperature is assumed 
along the way, rather than being derived. Properties of the theory are reported in the 



resuming tables in Appendix D 



5. Conclusions 



There are two viable perspectives on the results hereby derived. On the one hand, 
they might be seen as an exercise on abstract symmetry principles, later to be applied 
to theories that already enjoy diffeomorphism invariance, such as physics on curved 
spacetimes and General Relativity. We point out that relativistic diffusion is a broadly 
studied subject, and that it has some crucial peculiarities that do not allow direct 
application of our results. In this respect, although our Langevin equation is similar 
to Eq.(51) in Ref.jlH] (but for a term 2r>''P), their meaning is very different as to 
the role of time and the relationship between the diffusion matrix and the metric. Still, 
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covariant diffusion has an interest on its own. Graham introduced his own theory with 
these words: "Physical properties are independent of the coordinates used. Hence they 
must be formulated in terms of covariant quantities" ^U\. The will to "covariantize" 
(notions of) nonequilibrium thermodynamics was the ultimate drive for the present work 
as well. Recently Smerlak theorized that, along with propagation, it might be possible 
to "tailor" the diffusion of light in suitably engineered metamaterials |19]. Our proposal 
and his belong to the same paradigm, that spatially-dependent diffusion coefficients 
might act as an effective curvature of space(time). 

On the other hand, our work might be seen as a doorway to state-dependent 
diffusion even for mesoscopic systems, both at the level of interpretation and for physical 
modeling. By no means do we claim that the generally covariant Langevin equation is 
the correct equation for state-dependent diffusion, but it should be in all situations 
where the statistical properties of noise do not dependent on one preferred coordinate 
system. This demanding hypothesis has to be addressed independently, for example 
when performing microscopic derivations. A class of systems where curved Brownian 
motion is relevant is protein transport inside and across cell membranes under the effect 
of membrane curvature (see the introduction of Ref.[37j and references therein). In a 
parallel publication we also argue that it might be relevant for Brownian motion in 
non-uniform temperature environments [18]. The latter conjecture is already testable, 
as experimental works by Langon et al. [21] and by Volpe, Brettschneider et al. [22] 
(after a few issues have been resolved [23]) display great control over Brownian particles 
under the effect of spatial dependence of the diffusion coefficient, in this case due to 
the presence of walls. It might be feasible to design similar experiments where spatial 
dependence is due to temperature gradients or both to temperature gradients and spatial 
inhomogeneity, and put several of the models we described above at work. The generally 
covariant Langevin equation could also find application in Metropolis Langevin Monte 
Carlo sampling methods, according to the lines of research developed in Ref . |50] . 

The relationship between gauge invariance and detailed balance was a main theme 
in this article. It follows from our discussion that gauge invariance is necessary 
(but not sufficient) to write thermodynamic theories of state-dependent diffusion that 
enjoy equilibrium steady states, as different choices of gauge will otherwise lead to 
nonequilibrium driving. This was the case for Graham's theory. Surprisingly, instead, 
curvature does not affect the equilibrium character of steady states; it might be 
interesting to consider instead the thermodynamic effect of geometric torsion. 

Finally, we point out that most thermodynamic considerations in this work were 
carried out at the macroscopic level of the Fokker-Planck equation, and regarded only 
the equilibrium/nonequilibrium character of the steady state and the covariant definition 
of thermodynamic quantities. A great deal of covariant stochastic thermodynamics, in 
the spirit of Ref. [51] , has yet to be developed. 
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Appendix A. Exponential decay to equipartition 

We show that average energy has an exponential decay to its steady value, given by the 
equipartition theorem. We consider energy as defined in Eq.( l22l) — setting mA^ = 1 — 
and take its stochastic differential. The key technicality is to use Ito's Lemma 

dH = l^rfx" + ^dv'' + -^^^dv'^dvK (A.l) 

where we only retain contributions of order dt. Plugging into the generally covariant 
Langevin equation we obtain 

g LT g2 LT f)]-f / \ 

dH = —v^dt + ——-6'^ae1e]dt + — (--iv^dt - Tlv^v'dt + ^ eUW'\ . (A.2) 
Standard manipulations lead to 



E = -27 [H-'^j+ V2^g,,v^e\C. (A.3) 

Interestingly, the drift term only depends on the energy itself, so that we can average 
out the noisy term and obtain a closed equation for the average energy Et = (Ht), with 
solution 

It decays exponentially to its steady value 77,(7/27, which by the equipartition theorem 
is identified with nksT /2, yielding the Einstein relation. 



Appendix B. Overdamping limit of metric/diffusion theory 

We show that the generally covariant stochastic differential equation Eq. fl28|) (Brownian 
motion on a manifold) and the generally covariant Langevin equation Eq. ffTOj) are 
equivalent in the overdamping limit a —t- 00, at fixed temperature. More precisely, we 
work with the corresponding diffusion equations, from Eq. (!20|) to Eq. (l27I) . In fact, we 
prove the result for the more general theory (G) with independent metric and diffusion 
matrix discussed in SecHl and specialize the result at the very end. The treatment traces 
out Schuss's [52], §8.2]. To simplify, solutions to the perturbative chain of equations 
that will emerge are given, as they are known from previous literature. A logically close 
step-by-step derivation can be operated by means of expansions in terms of Hermite 
polynomials, as is done in Ref . [25] . 

Set both /3 = 1 and [2TrhY / mc^ = 1. Define the adimensional energy 



(B.l) 
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Let e = a~^. We rescale time according to t = e~^s in the covariant Kramers equation 
Eq.f l20|) . and expand to order e both the generator and the probabihty density: 

(^e^L2 + eLi + Lq) (Pq + ePi + e'Ps + • • • ) = 0, (B.2) 

where 

(9 / dP\ 

LoP = 9^'g^^[h,v^P+g^,), (B.3a) 

dP B 

L^P= -'^''o^a + g^ai^lv'v^P)^ (B.3b) 
dP 

L,P = - (B.3c) 

We obtain the following expansion: 

= LoPo, (B.4a) 

= LoPi + LiPo, (B.4b) 

= L0P2 + LiPi + L2P0, (B.4c) 



A solution of the first equation is 

Po(x, s) = qo{x, s) exp -0(a;, w), (B.5) 

where qQ{x, s) is an undetermined function. The spatial density is found by performing 
the integration over velocities: 

Pq[x,s)= I dv Po{x,v,s) = — 7==^- (B.6) 
J \/h[x) 

We also define the conditional probability 

7r{v x) = — = exp-0 (B.7) 

Po{x, s) N 



which is a multivariate normal distribution. Plugging Pq into Eq.( 1B.4b| . we obtain 



A solution is given by 

Pi = (gi-^;"^,)Po, (B.9) 
where gi = qi{x, s) is again an undetermined function and 6a is easily seen to be 

9a = gabh''g^Jog^. (B.IO) 
Substituting into Eq. flB.4cp and integrating with respect to velocities gives 
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where it is assumed that all terms decay sufficiently fast for high velocities, so that 

Eq. dRTT]) then yields 

„ d 



dpo 
ds 



dx°- 



b d , go , D 



dv. 



(B.13) 



The ffist integrand is odd and its contribution vanishes. In view of Eq. flB.6p . we obtain 



dpo d 



ds dx"^ 
We now perform the integration 



dc 



d 



Po 



v"'v Po{x, V, s) dv = po{x, s) / v'^v 7r{v\x)dv = po{x, s)h°' (x) 



(B.14) 



(B.15) 



(B.16) 



In particular, when hab = Qab we obtain Eq.( l27l) . It follows that in the overdamping limit 
the generally covariant Langevin equation and brownian motion on a manifold furnish 
equivalent descriptions, up to ffist order in e. 



where we recognized the covariance matrix of -7r(-,x). Finally we obtain 

ds dx^ \ dx^ ^/h 



Appendix C. Overdamping limit of diffusion with state-dependent 
temperature and viscosity 

We perform the overdamping limit of theory (F), with state-dependent temperature and 
viscosity. The particular case with state-dependent temperature is obtained by setting 
l^ab _ gab^ rjj^^ main difference from the above overdamping limit is that the term with 
Christoffel coefficients disappears, resulting in a minor complication. All goes as before, 
but for the definition of the operator 



LiP = -v" 



dP 

dx°- 



(C.l) 



The solution to the ffist equation in the pertubative chain is the same. As above, 
plugging Po into Eq.( IB.4b1) . we obtain 



gab 9 f l^^^^cp^ ^ ^ 

dv°- \ dv^ 
We look for solutions in the form 



d 



v''—logqo-v%'v''natc]Po- 



Pl={qi + ^V^'v^'J^abc - V^Oa ) Po, 



(C.2) 



(C.3) 
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where qi = qi{x, s) is again a scalar function, 6a is a vector field and T^abc a tensor. After 
some work we obtain that T,abc = Qad.h'^'^^ebc and 



Oa = gabh'" ( ^ log go - ^h'^^n^ea - ^/^'^^^ade 



../.-(^elog,0-^^.logv/)I+^/..../.-) 

(C.4) 



make Eq. (lC.2p satisfied, where we used the identities 

h'^nbca = da log Vh, h'^nabc = -hacObh'^ " Oa log Vh. (C.5) 

Substituting into Eq. (lB.4cP and integrating with respect to dv we obtain 

% [ dv (-v''v% + \v%%%^gbfhf'n,,,] Po ^ 



dpo 
ds 



where again the term with qi drops for it gives an odd integrand. Given Eq. (lB.15p . the 
first current contribution reads 

Ji = Poh^'gbc (^-d' log Po-^d'^ log ^/h - ^d.h"''^ . 

where we raised indices with h"'''. By the Isserli-Wick theorem 

j dvv'^v^v^v'^Po = Po {h^^h"^ + h'^^h!"^ + h'^^h^") . 

Letting G"^ = h'^^g^^h'^^ 



{3I+32) (C.6) 



(C.7) 



(C.8) 



32 



Po 
3 L 



h^^gbci-ddh'^" - d" log \fh) + 2h''^G'"^a 



cdb 



for a total Fokker-Planck current: 



3 



d"Po+ log + 9^/1'^^^ 



Po 



(C.9) 



(C.IO) 



In particular when /i"'' = g"''^ we obtain the Ito generator dsPo = —dbda [g"'^Po) , while in 
the 1-dimensional case we find: 

dpo 



ds 



9_ 



-dpo + ^{2 + ^]d\ogh 



(C.ll) 



Appendix D. Comparison of models 



The following diagram summerizes the interdependence of the several Langevin 
equations that we considered. The symbol -w means "is replaced by" . 



(D.l) 




Resume tables can be found in the next page. 



Table Dl. State-dependent diffusion theories in second-order formalism. 





Model 


Refs. 


Langevin equation 




Einstein relation 


Steady state 


A 


Generally Covariant 


ivi, [361 EZ] 


x"" + T'^^x''x^ = — -|- V 




7 = /3o- 


fifexp -^gabv'^v'' 


B 


Covariant (Graham) 


[101 [m 






7 = /3(T 


? 


C 


State Dep. Temperature 


[Ml ESI EHl EH EH] 






7 = /3cr 


7 


D 


State Dep. Viscosity 


[m EQl sal E] 


x« = -^ix^ + ^elC 




7," = Pag'^^b 


exp — |(5afef 


E 


Pure posticipatory SDE 


[51 Sec.V] 










F 


St. Dep. Vise, and Temp. 








7b" = Pag'^'h.b 


? 


G 


Metric /Diffusion 


ivi 


x"" + n'^^x^x'' = -7^ x'l + 




^b = Pag'^'Kb 


hexp —^habV^v^ 



ilL are the Christoffel coefficients of h"-^ 



Table D2. State-dependent diffusion tlicories in first order formalism. Symmetries and thermodynamic properties. 





Smoluchovski eq. 


FP current 




ID S.S." 


Gauge inv. 


Covariance 


D.B.^ 


A 




-g"^ {dbP-pdblog^) 




Yes 


Yes 


Yes 


B 




-g^^dbP - {6'^efdbe'') p 


V9 


No 


Yes 


No 


C 


x" = 


-dbig^'p) 




9 


Yes 


No 


No 


D 


i« = dbg''^ + ^eK' 


-g"'''dbP 




const. 


Yes 


No 


Yes 


E 


x" = 2 5'^e%e] + V2e^C 


-g''^dbP + 5'^ {epbef - efdbc'^) p 


const. 


No 


No 


No 


F 


(omitted) —h"-^gbc 


d^p + (id' log ^^h + ddh'^'') p 


+ lh'^^G'''n,dbP ^ 


2 r gdh 

h3e-' 3^ 


Yes 


No 


No 


G 


^ -gbc^a^ - df, log ^yh/g + ^2 e'^C 'Q"^ (Obp-p db log Vh^ 




Yes 


Yes 


Yes 



° Steady state in the 1-dimensional case. 

^ Detailed Balance, independently of the specific choice of metric, coodinates and gauge. Specific models might obey detailed balance. 

t Qab ^ ^acg^^f^db^ 
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